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We study the problem of master-slave synchronization and control of totalistic cellular automata 
(CA) by putting a fraction of sites of the slave equal to those of the master and finding the distance 
between both as a function of this fraction. We present three control strategies that exploit local 
information about the CA, mainly, the number of nonzero Boolean derivatives. When no local 
information is used, we speak of synchronization. We find the critical properties of control and 
discuss the best control strategy compared with synchronization. 



I. INTRODUCTION 

There is a certain interest in modelling extended sys- 
tems using discrete units, in many cases Boolean ones. 
There are two main areas for which discrete modelling is 
appropriate: systems whose dynamics are well approxi- 
mated by a discrete process, and those for which discrete 
processes constitute a simplified computational descrip- 
tion. Examples of systems of the first kind are genetic 
networks [1], some formalization of neural networks [2], 
DNA replication and translation [3] , and VLSI digital cir- 
cuits. These discrete models are stable and predictable 
in spite of a noisy environment. In order to overcome 
noisy fluctuations affecting signals of electronic devices, 
a large negative feedback is added, decreasing the over- 
all efficiency but making the system more reliable. The 
extreme case of this approach is that of working at the 
saturation level of devices (e.g., transistors), leading to a 
digital representation of the problem. In biology, DNA 
replication and translation is essentially a digital pro- 
cess, and in many aspects gene control is also a kind of 
deterministic device, even in the presence of a fluctuating 
environment. 

On the other hand, discrete modelling and in particular 
cellular automata are widely used as effective computa- 
tional tools for simulating complex systems [4] . The idea 
is that of looking for the emergence of a complex dynam- 
ics out of an ensemble of simple components. Some exam- 
ples are opinion formation, epidemics, and traffic [5, 6]. 

A discrete system is formed by separate units which 
can be nodes or sites. The state of each unit is taken 
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from a discrete set and evolves in discrete time steps. 
Cellular automata (CA) [7] are the typical mathemati- 
cal examples of such systems, even though one may be 
also interested in non-homogeneous contact networks and 
non-parallel dynamics [8]. 

We aim to study the dynamical properties of discrete 
dynamical systems, and in particular how these charac- 
teristics may be exploited for their control. In the fol- 
lowing we shall consider simple one-dimensional CA, but 
most of our considerations can be extended to more com- 
plex systems. 

The state of a discrete dynamical system is given by a 
discrete set of numbers, and a trajectory is simply given 
by a list of such states. Since the total number of possible 
states is finite, and due to the deterministic character of 
the system, after an eventual transient the trajectory has 
to enter a cyclic behavior (including fixed points). How- 
ever, the size of this periodic attractor and the number 
of different attractors may scale with the size of the sys- 
tem. These numbers may grow exponentially with the 
system size, and quickly become intractable for practi- 
cal purposes. We shall denote these as discrete extended 
systems. 

Discrete extended systems are not chaotic in the usual 
sense. They have been selected or built with the explicit 
goal of being insensitive to small perturbations. How- 
ever, they react to finite stimuli (change of state of some 
element), as in the experiments of gene knockout or si- 
lencing [9]. In this case, the reaction may be localized, 
or may extend to the whole system. In the latter case, 
the future evolution of the system is practically unpre- 
dictable. This is a situation analogous to usual chaotic 
systems, and the term "stable chaos" has been coined for 
such systems [10]. 

Although classical perturbation theory cannot be ap- 
plied to these systems, it is possible to define the discrete 
analogous of derivatives (Boolean derivatives [11, 12]), as 
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illustrated in Sec. II. By using this concept, it is then 
possible to construct a Jacobian matrix and to define the 
largest Lyapunov exponent, which discriminates well be- 
tween systems that are "stable chaotic" and others whose 
evolution is more predictable [13]. 

A different, but complementary, approach to the study 
of the chaotic properties of a system is based on master- 
slave synchronization [14]. The idea is that of driving 
the slave with part of the signal from the master. It is 
possible to show that, for simple maps or for a uniform 
driving in extended continuous systems, there is a rela- 
tionship between the synchronization threshold and the 
largest Lyapunov exponent. 

A modified version of this synchronization procedure 
can be applied to discrete systems. The idea is that of 
synchronizing them using an all-or-nothing procedure; a 
fraction of sites are chosen at random and the state in 
the slave system at those sites, is set equal to that of of 
the corresponding sites in the master system (pinching 
synchronization) [13]. Applied to cellular automata, the 
pinching synchronization threshold shows a remarkable 
correlation with the previously defined largest Lyapunov 
exponent defined using Boolean derivatives. 

The pinching synchronization procedure defines a di- 
rected percolation process among defects, that is, the 
sites that differ between the master and the slave. We 
say that the two replicas are not synchronized if in the 
long time limit there is some defect surviving (a more 
precise definition is given in the following section). 

The synchronized state is absorbing, once it is reached, 
it cannot be abandoned. However, in the absence of a 
synchronization effort, the synchronized state is unsta- 
ble for "chaotic" systems, since a finite perturbation will 
in general propagate to the rest of the system (as said, 
this is indeed a possible definition of a chaotic discrete 
system). The synchronization efforts (pinching synchro- 
nization) has the goal of making the synchronized state 
stable, and the synchronization transition corresponds to 
the marginal stability of this state. 

By considering now only the points to which the pinch- 
ing synchronization is not applied, one has a percolation 
backbone over which the defects can survive. Clearly, 
there is an upper value of the pinching synchronization 
probability over which the percolation backbone disap- 
pears and and synchronization takes place regardless of 
the rule. Differently from highly chaotic maps [15], all 
cellular automata synchronize well below such limit, in- 
dicating that the self-annihilation of defects plays a fun- 
damental role in the synchronization transition. Indeed, 
the disappearance of defects is due both to the pinching 
procedure (which can be random and uncorrelated) and 
to the self-annihilation process, which in general is cor- 
related. This correlation has a deep impact on the criti- 
cality of the phase transition, as we show in Section III. 

We finally arrive to the main point of this paper. In 
control theory a controller is used to make a dynamical 
system behave in a specified way. The problem of control 
of a dynamical system may be split into two parts: the 



driving of the system to a target area in phase space, and 
the stabilization of a given trajectory originating from 
this area. Chaotic systems are ideal targets for control, 
since their sensitivity to small changes may be exploited 
to drive them to the target area [16]. After that driving 
phase, chaos may be suppressed in order to make them 
follow, for instance, a desired periodic orbit [17]. 

We shall concentrate on the second part, the stabiliza- 
tion of a given trajectory. Let us denote with the term 
"natural" a trajectory that is generated by the dynamics 
of the system, and with the term "artificial" all the other 
trajectories. In other words, if we initialize the given sys- 
tem in one of the states of a natural trajectory, then the 
system will follow the trajectory without any additional 
external control. On the contrary, "artificial" trajecto- 
ries need a continuous control to force the system under 
study to follow them. 

The problem of driving a chaotic system on a natural 
trajectory may be seen as the problem of synchronizing 
a "slave" with a "master" system that happens to fol- 
low the desired trajectory. However, while in the studies 
about synchronization one exerts little attention to the 
optimization of coupling, when formulated as a control 
problem this becomes a crucial issue. 

In synchronization, the effort is applied at every site. 
In control problems, one wants to exploit some knowledge 
about the system. It is therefore analogous to a synchro- 
nization problem of two different systems with a "tar- 
geted" force, that tries to "kill" the defects in an efficient 
way. We show how the concept of Boolean derivatives 
can be exploited to achieve this goal [18]. 

The result is however rather surprising: the control 
efforts should concentrate on the regions giving birth to 
a small number of defects, while chaos can be exploited 
to "self-synchronize" the systems. 

In Section II we present definitions of totalistic cellu- 
lar automata, of Boolean derivatives, and synchroniza- 
tion and control of CA. Results of synchronization and 
three control strategies are discussed in Section III and 
in Section IV we present some final comments. 



II. CHAOS, SYNCHRONIZATION AND 
CONTROL OF CELLULAR AUTOMATA 

In this Section we recall the definition of totalistic 
Boolean cellular automata of range R, of Boolean deriva- 
tives, and present the master-slave synchronization and 
control of cellular automata [12, 13, 19]. 

A Boolean cellular automaton of range R is a map 
/ : B N -> B N with B = {0, 1}, N large, and / = 
(/o,...,/jv_i) with fi : B N B, i = 0,...,N - 1. 
A state of the CA is x = (x , ■ ■ ■ ,xn-\) with Xi G B, 
i = 0, . . . , N — 1. The map / defines a flow on B N , 

x(t+l) = f(x(t)) (1) 

with t = 0, 1, In what follows /, = /, Vi. 
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Let Mi be a neighborhood of site i with R sites and 



with < Si < R. For totalistic CA of range R, 

Xi(t + 1) = f(Si(t)), 



(2) 



(3) 



with / : {0, . . . , R} -> B. There are 2 R+1 different CA of 
range R, each one defined by a (R + l)-tuple (yo, ■ ■ ■ ,Ur) 
such that yj — 1 if the outcome of the CA is one when 
Si = j and yj = otherwise. In what folows, we use 
Vichniac's notation RTC, with R the range, T is for to- 
talistic, and C the number in base 10 of (yo, ■ ■ ■ , Vr) [11]. 

The control mechanisms wc consider depend on the 
knowledge of the local expanding properties of the CA 
given by the Boolean derivatives [12] defined by 

T _ dfi(x) 

=fi{x , ...,Xj®l,.. . ,X N -i)® 
fi(x , ...,Xj,.. .,X N -i) 

{1 fi changes when Xj changes, 
fi docs not change when Xj changes, 

with © the logical exclusive disjunction or what is the 
same, the sum modulo 2. 

Higher order Boolean derivatives can be defined in a 
similar way and we can write any function / of range R 
in the equivalent of a Taylor or Maclaurin expansion that 
is exact if the expansion goes up to i?-th order deriva- 
tives [12]. The exact MacLaurin expansion of a CA with 
R = 2 is 



f(x,y) =/(0,0) 



df 
dy 



dx 
■y® 



x® 



(0,0) 

d 2 .f 



(0,0) 



dxdy 



(0,0) 



(4) 

■x-y (5) 



with • the logical conjunction. If the expansion of / con- 
tains only the first order Boolean derivatives, we say 
that / is linear. As an example, let R = 3, A/i = 
{xi-i, Xi, Xi+i}, and /(sj) = 1 if Sj is odd and zero oth- 
erwise. This is the 3T10 CA shown in Table I which is 
linear since 

f(Xi-i,Xi,X i+ i) = Xi-i ®Xi ®x i+1 . 

In Wolfram's notation, this is rule 150 [20]. Another 
example with R = 3 is the 3T6 CA, rule 126 in Wolfram's 
notation, defined in such a way that /(sj) = 1 if Si = 1, 2 
and zero otherwise shown in Table II. It is nonlinear 
since 

f(Xi-i,Xi, X i+ i) =Xi-i ® Xi® x l+ i® 

Xi—iXi © Xi — \Xi+\ © XiXi-\-\. 



TABLE I. Truth table and Boolean derivatives of the 3T10 
CA. 
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TABLE II. Truth table and Boolean derivatives of the 3T6 
CA. 
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We now discuss a master-slave synchronization and 
control mechanism of cellular automata where the con- 
figuration x represents the master and y the slave. Let 
Pi, i = 0, . . . , N — 1 be a Boolean variable that is 1 with 
probability iTi and otherwise. That is, 



Pi = 



1 with probability 7r^ 
with probability 1 — Hi 



Then 



Xi (t + l) = f(x(t)), 

Vl {t + 1) = (1 - Pi {t)) ■ f(y(t)) © Pi (t) ■ f(x(t)), 
= {f{x{t)) if Pi (t) = l 

\f(y(t)) ifPi(t) = o, 
^ + 1 ) = w(* + i)eii(t + i). 

In the last Eq. hi = 1 if yi is different from Xi and is zero 
otherwise. We use the shorthand h = y © x. Wherever 
Pi = 1, y% — x i and we say that y and x are "pinched" at 
site i. We define the control parameter k and the order 
parameter h by 



k 



N 



(6) 



If k = 0, y evolves independently of x and h ^ 0, even 
though the self-annihilation of defects is still present due 



FIG. 1. Time evolution of the difference h for (a) CA 3T10, 
(b) CA 5T42 with N = 200 during 200 time steps. Time runs 
from top to bottom. 
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FIG. 2. The scaling of the asymptotic distance h as a function 
of the synchronization effort k for rule 3T10 in case and N = 
20, 000. The curves follow the scaling law Eq. (8). The dotted 
curve corresponds to ft = 0.27 (k c ~ 0.321), compatible with 
directed percolation and the dashed line corresponds to /3 = 1 
(k c — 0.3275) compatible with a mean field behavior. In the 
inset, the scaling behavior is reported on a log-log scale. 

to the Boolean character of the system. If k = 1, h = 0, 
both states are pinched together. There is a threshold k c 
above which h = in the long time limit. See Fig. 1 for 
a visual illustration of the synchronization process. If no 
knowledge is assumed about the local behavior of / that 
defines the CA, we speak of pinching synchronization [13] 
and take 71"$ = p for all i. In the case where some knowl- 
edge about / is used to determine 7Tj (that can change in 
time), we speak of pinching control. 

In the following Sec. we show results for synchroniza- 
tion and three types of control based on the local expand- 
ing properties of the map /. We define Ji as the sum of 
the R Boolean derivatives that are affected by site i 

Ji — ^ ] Jj.ii CO 

with < Ji < R (notice the indices). The cases we 
consider are 

0. iTi(t) — p with < p < 1 independently of i and t. 

1. 7r,(t) is proportional to Ji(t). 



o 1 1 1 1 1 1 1 

0.05 0.1 0.15 0.2 0.25 0.3 

k c (R) - k 

FIG. 3. Synchronization (control type 0) of some linear cel- 
lular automata, from top to bottom: 3T10, 5T42, 7T170, 
9T682, rescaled to make the critical point coincide with the 
origin of axis. In these simulations N = 4, 000. 
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FIG. 4. The relation between k c , and the range 7? for linear 
cellular automata, control type 0. The solid line marks a 
the best fit of Eq. (??) with k^, = 0.825, 7 = 0.282, and 
a = 0.828. 

2. If Ji(t) = 0, TTi(t) — 0, otherwise 7Tj(t) is propor- 
tional to R — Ji(t) + 1. 

3. If Ji(t) = 0, TTi(t) — 0, otherwise 7Tj(t) is propor- 
tional to (R - Ji(t) + l) 2 . 

Case corresponds to synchronization, no knowledge on 
the dynamics of the CA is used. The other cases, use the 
derivatives in one way or the other. Since < Ji < R, 
case 1 will favor that states x and y be pinched at site i 
at time t when Ji(t) is large. On the other hand, cases 2 
and 3 favor pinching at site i at time t whenever Ji(t) is 
small (case 3 more than case 2). 



III. RESULTS 

We begin by presenting the synchronization and con- 
trol of linear CA. All types of control coincide with the 
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FIG. 5. Critical behavior of nonlinear CA 7T10 for control 
types 0, . . . , 3. In all cases Eq. (8) fits the data with ho ~ 0.46 
and different values of k c and /3 (continuous curves). Control 
type and 1 agree numerically with j3 ~ 0.36, control type 2 
with /3 = 0.32 control type 3 with = 0.28. 

synchronization case for linear CA, since Ji = R for 
every site i. 

Near k c , h follows the scaling law 

with h the average number of active sites for the unper- 
turbed rule (k = 0). 

In Fig. 2 we show the critical behavior of h for CA 
3T10 that has a crossover between the directed perco- 
lation (DP, (3 ~ 0.27) and the mean field (MF, (3 = 1) 
universality classes. As we show in Fig. 3, other linear 
rules exhibit the same type of crossover from the MF to 
the DP universality class. As R grows, the MF character 
of the transition occupies a smaller region. This behavior 
is rather unexpected, since the Grassberger- Jensen con- 
jecture [21] should hold in this case. The interactions are 
short-range and, since the visual behavior is irregular, 
correlations are expected to be short-range. By applying 
a transformation y = 1 © x for odd neighbourhoods one 
can show that for a linear CA, ho = 1/2. 

In Fig. 4 we show the synchronization threshold k c of 
linear totalistic CAs as a function of the range R. The 
numerical data are well fitted by 

fc c = fc 00 (l-cxp(- 7 (i?-l)«)). (9) 

We assumed that for R = 1 the critical value of k c is 
zero. In this case each site only has one neighbor, and 
therefore any infinitesimal control is able to synchronize 
the replicas in the long time limit. When R — > oo, k^ ~ 
0.825. 

For the nonlinear CA 7T10, the scaling behavior near 
the transition is more similar to a DP character, for all 
types of control, as shown in Fig. 5, even though the 
measured exponents differ numerically from the DP ex- 
ponent. 



7T6 7T30 3T6 




6T60 5T44 5T46 




5T45 6T28 7T60 




FIG. 6. Phase diagram of several nonlinear CA with N = 
4,000. 

For linear CA, all control cases are equivalent. This 
is not the case for nonlinear CA. Although we have not 
performed an exhaustive search, some nonlinear chaotic 
CA show the general pattern k c < kc < k c ^ < k^ 

(n) 

with kc the critical value of the control parameter k in 
case n, n = 0, . . . , 3. Examples are CA 7T6, 7T30, and 
3T6 whose critical behavior is shown in the top three 
graphs of Fig. 5. For the CA of the second line of Fig. 5 
control 3 shows the smallest k c . In the last line of Fig. 5 
we show three nonlinear CA for for which the four types 
of control are similar. 

In a large sample of nonlinear CA studied, that in- 
cludes those of the previous Fig. , the application of con- 
trol 3 gives the smallest value of k c (which may coincide 
with the other types of control) . A possible explanation 
is that, due to the irregular behavior of the CA and the 
presence of an absorbing state (the synchronized or con- 
trolled state), the system will approach this absorbing 
state in small neighborhhoods that can grow at most lin- 
early in time. Thus, it is sufficient to break down the 
system into small portions, and for doing this it is more 
efficient to synchronize sites with a small (but non-zero) 
number of derivatives, which are presumably on the bor- 
der of the non-synchronized patches, than synchronizing 
other sites in the bulk. Control 2 and even more control 
3, preferably select sites with few non-zero derivatives, 
and are therefore more efficient. On the other hand, 
it is almost useless to synchronize sites in the bulk of 
the patches, with many derivatives, since they do not 
contribute much to the spreading of the desynchronized 
state. 

This is not the whole story. As we show in Fig. 6, there 
are CA for which control type 3 is not the best one, even 
if the differences are minimal. At present, we have no 
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explanation for this behavior. 



IV. CONCLUSIONS 

We have studied the problem of master-slave synchro- 
nization and control of totalistic Boolean cellular au- 
tomata via a pinching (all-or-none) mechanism. We have 
investigated four different cases: pinching synchroniza- 
tion, 1 control proportional to the number of Boolean 
derivatives at a given site, and 2 and 3 controls for which 
the probability of applying the synchronization proce- 
dure diminishes with the number of Boolean derivatives. 
The number of nonzero local derivatives is an indicator 
of the expected number of defects. 

For linear rules, the four types of control are equivalent 
and the control transition shows a crossover from a di- 



rected percolation universality class to a mean field one. 
More research is needed to decide this, since It is well 
known that the estimation of the exponent j3 requires 
sophisticated techniques which are left for future study. 

Nonlinear CA show a counter intuitive behavior. For 
some of them, control 3 is the most efficient followed by 
controls 2, 0, and 1, although there are CA for which 
this is not true. A more exhaustive research is needed 
in order to find some way to group CA according to the 
order in which they can be controlled. 
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